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Abstract By applying hybrid functions of general block-pulse functions and the second 
Chebyshev polynomials, integrodifferential systems are converted into a system of algebraic 
equations. The approximate solutions of integrodifferential systems are derived. The nu- 
merical examples illustrate that the algorithms are valid. 
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1. Introduction 

Integrodifferential systems are a class of systems of importance [T]. Since it is difficult 
to obtain the analytic solutions of integrodifferential systems, numerical methods to obtain 
approximate solutions are interested. Many orthogonal functions or polynomials, such as 
block-pulse functions [21 [3] , Walsh functions [5] , Fourier series [5] , Legendre polynomials 
[6], Chebyshev polynomials [7j and Laguerre polynomials [8], were used to derive solutions 
of some systems. In recent years different kinds of hybrid functions [TU] and Legendre 
wavelets were used. In this article we apply the hybrid functions consisting of general 
block-pulse functions and the second Chebyshev polynomials to converting integrodifferen- 
tial systems into a system of algebraic equations. The hybrid solutions of integrodifferential 
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systems are obtained. We present the general operational matrices. The numerical examples 
show that this method is convenient for application. 



2. Preliminaries 

2.1. Definitions 

A set of block-pulse function bk{t), k = 1,2, . . . , K, on the interval [to,tf) are defined as 

i 0, otherwise, 

where tx = tf and [tk-i,tk) C [to,tf), k = 1,2, . . . ,K. 

The m order second Chebyshev polynomials in the interval [—1,1] are defined by the 
following: 

Sm{t) = sin((m + l)arccos(t))(l - t^)"^, (1) 
which are also given by the following recursive formula 

r So{t) = l, 5i(t) = 2t, 

I S,n+i{t) = 2tSm{t) - Sm-l{t), m = 1,2, 

The hybrid functions hkm{i), k = 1,2, . . . , K; m = 0, 1, . . . , M — 1, on the interval 
[to,tf) are defined as 

hkmit) = bk{t)Sm{dkH2t - tk-l - tk)), (3) 
where dk = tk - tk-i, k = l,2,...,K. 

2.2. The operational matrix 
Let 

Hk{t) = [hko{t), hk,M-i{t)V, H{t) = [Hl{t), Hj,{t)Y. 

It is not difficult to verify that 
't 



[ H{t)dtc^PH{t), 

J to 



(4) 



where 



K-lK-k 



fe=l i=l A;=l 



^ k=l k=2 ^ k=3 



E^J^^ is the mxm matrix with 1 at its entry (i, j) and zeros elsewhere, denotes Kronecker 
product and [^^y^] is the greatest integer part of ^^^^f^- 

2.3. Function approximation 

An Z-dimensional vector function f{t) on the interval [to,tf) is expressed as 

K M-l 

fit) ^Y.Y1 fkmhkmit), (5) 
fe=l m=0 

where 



fkm = - f f(P'~\dkt + tk-1 + tk))Smmi - t^f^dt. 

Rewrite f{t) as 

K 

f{t)c^^FkHkit) =FH{t), 
k=l 

where 

Fk = [fkO, ■ ■ ■■> fk,M-l\, F = [-^1, • • • , Fk\- 

For corresponding Fk and F we denote 

Fk = [/feO' • • • ' fk,M-iY ^ F = [Fl, Fj^Y, 
where r is the transpose. 

2.4. Expression of the product of a matrix function and a vector function 



(6) 
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Let a matrix function M{t) be appropriate to a vector function f{t). We express M{t) 
and f{t), respectively, as 

K M-l K M-l 

M{t) Mkmhkmit), fit) fkmhkmit). 

k=l m=0 k=l m=0 

Then 

K M-l M-l 

M{t)f{t) Mkifkjhki{t)hkjit). 

k=l i=0 j=0 

Prom 

M-l 



hki{t)hkj{t) ~ ^ %m^km{t), 



where 



= 1 r s,{t)Sj{t)Smmi - t't^dt, (7) 

J -1 



we have 

K M-l y M-l M-l 



M{t)m ^ E E ( E E 4^lMk^fk^hkm{t) 

k=l m=0 ^ i=0 j=0 ' 
K M-l K 
= E E fkmhkmit) =^MkHk{t), 



k=l m=0 k=l 

where 

M-l M-l 



(8) 



Mkm "EE '^kL^^kifkj — Mk^Fj;, Fk — [/fco, • • • , /fc,M-l]^' ^k — Mj^Fj., 



Mkm 

Therefore 



M-l M-l 
^ i=0 1=0 



K M-l 



N{i,t)f{t) C^Y^Yl ^km{t)hm{t), (9) 



k=l m=0 



where 



Nkm{t) 



i=0 



1=0 



Nki{t) = - [ N{i,2-\dkt + tk-i + tk))Si{t){l-t^)'2dt. 



Let 



j-tf _ K M-l 

(t) = / N{t,t}f{t)dt^Yl Yl Wkmhkmit), 
J to 1 1 ™_n 



k=l m=0 



K M-l 



Nkrn{t)c.Y,Y.^kydt), 



where 



(ji) 



km 



k=l m=Q 



M-l 



/ . "'km ^^ki ' • • • ' 



A/-1 



{i,M-l)Uil) 
km ki 



j=0 



1=0 



iV, 



ki 



Nkii2-\djt + tj^i + tj))Siit){l - t'y^dt. 



-1 



By Eq.Q we have 

K M-l K M-l 



w 



j=l 1=0 k=l m=0 



hji{t)dt]Fkhji{t). 

tfc-i / 



So 



[^] if K 



fc=i i=i j=i 



where 



AT 



j,2k-l 



N 



N 



{i,M-^V 



j,2k-l ' • • • ' ^^j,2k-l 



3. Analysis of integrodifferential systems 



Consider the following integrodifferential system described by 

(13) 



x{t) = A{t)x{t) + / N{t,t)x{t)dt + B{t)u{t), te[to,tf], 



L x{to) = Xq, 

where x{t) is an n-dimensional vector function and u{t) an r-dimensional vector function. 
A{t), B{t) and N{t,t) are the matrix functions of appropriate dimensions. Express A{t), 
x(t), J^f N(t,t)x{t)dt, B{t) and u(t) respectively, as Eqs.® and i^. By Eqs.® and © 
we have 

K M-l K M-l 

A{t)x{t) ^^Yl ^kmhkmit), B{t)u{t) ^ Yl ^kmhkm{t), 

k=l m,=0 k=l m=0 

tf K M-l 



N{t,t)f{t)dtc^Yl Yl ^kn,hkm{t). 



*0 k=l m=0 

Integrating Eq. (fT3]) from to to t and combining Eqs.®, Q and (fT2ll we obtain 
[Xi, . . . ,XK]H{t) — [Xoi, . . . ,XoK]H{t) 

= f\[Ai, ...,Ak] + [Wu..., Wk] + [Bu BK]}H{t)di 

J to 

where 

Xk = [xko, ■ ■ ■ ,Xk,M-i], Xok = [xo,0, ... ,0], k = 1,2,. . . ,K. 

Thus 

[Xi, ...,Xk]- [Xoi, ...,Xok] = {[Ai, ...,Ak] + [Wi, . . . ,Wk] + [Bi, . . .,Bk]}P. 
Using Kronecker product we rewrite the above equation as 

x-xo = iP-^ inKiM, hr + [wi, wkv + [ir, . . . , hv}, 

where 

^ = 1^1 1 • • • ) ^kVi = [XkO^ ■ ■ ■ ) ^k,M-lV^ = [^fcO' ■ ■ ■ ' ^k,M-lV' 

Xo = [Iji, • • • , X^kV, Xok = [xl, 0^ . . . , OT, k = l,2,...,K. 
Wk and Bk have the similar meaning as Ak- So 

X-Xo = {P^® . . . , {AkXkYY + W+ [{BiUiY, {BkUkYY} 
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K [^] K K 

k=l k=l i=l j=l 



fc=i ) 
+E(^, 



A:=l 



Therefore 



X = ru + n, 



where 



K 



r = [l 



MKn 



k=l 



K 



[^] X if 



^ = E(4. ^^^)+ E EE^^ 

A:=l fc=l i=l j=l 



(14) 



4. Numerical examples 



4.1. Example 1 



Consider the integrodifferential system described by 



x{t) 



t' + l -t 
1 



x{t) + 



t 3 
3t^ 



x{t)dt + 



2e - e 



u{t), 



x(o) = [o,o]^ 



where x{t) is a 2-dimensional state function and u{t) a 1-dimensional control function. 

Let u{t) = 1. Then by Eg. (1141) and choosing M = 4 and K = 3, we have the 
hybrid solutions x{t) = [xi(t) , X2{t)Y : 

Xiit) = ^/lio(t) + ^/lll(t) + m^uit) + ^/i20(t) + n/i2l(t) + Ti4/l22(t) 
+ Tii^30(t) + i/l3l(t) + Ti4^32(t), 



which can be checked to be the same as the exact solutions x{t) = [t'^,t^Y. 
4.2. Example 2 



Consider the integrodifferential system described by 

• 1 



x{t) 



1 t 
t t' + l 



x{t) + 



+ 



3t e-* - P 

3e-i - 5 - 3t 
2e-^ -7 -t-?,f 



x{t)dt 
u{t), x(0) = [l,3]^ 



where x{t) is a 2-dimensional state function and u{t) a 1-dimensional control function. 

Let u{t) = e~*. Then by Eg. (1141) and choosing M = 5 and i^T = 3, we have the 
hybrid solutions x{t) = [xi(t) , X2{t)Y : 



+ 



465 : 



764 ^2o(^) - 4007^21 (^) 

^30 (^) — ^)^3l(0 
613 



22756 ^^^^ 36694 "i'^v v ' 587104 
121 A /+^ 5 u U\ I 1 



1213 
2877 



h22it) - 



+ 

= ifM^lo(i) - ^hu{t) + ^^hi2{t) - 36|94^13(^) + 39^^14(^) 



31092 
187 



85351^23(1^) + 819370 
9 



/i24(t) 



+ ^h2Q{t) - ^h2l{t) + ^^h22it) - 7^/i23(i) + 54^47^24^ 



1744 L _ 449 L /^N I 47 l _ 

1333"'30l''J 4123"'31l'^/' "I" 10364 ""32 1*-; 63529 '"■^•^^"^ ' 381175' 



) '''33 

W + 381175 ^34 (t). 



Table 1 and Table 2 imply that the hybrid solutions converges the exact solutions 



x{t) 



r*,3e-T 



rapidly. 



6. Conclusion 

Using the excellent properties of operational matrices of the hybrid function of 
general block-pulse functions and the second Chebyshev polynomials, the general al- 
gorithms for integrodifferential systems are derived. The illustrative examples demon- 
strate that this technique is convenient for application. 
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Table 1: Analytic and hybrid values of xi{t). 



t 


Analytic 


Hybrid (Ar=4,M=5) 


Hybrid (Af=4,M=7) 


Hybrid (iV=4,M=9) 


0.1 


0.90483741803596 


0.90483741135846 


0.90483741803662 


0.90483741803596 


0.2 


0.81873075307798 


0.81873074696139 


0.81873075307870 


0.81873075307798 


0.3 


0.74081822068172 


0.74081822463206 


0.74081822068111 


0.74081822068172 


0.4 


0.67032004603564 


0.67032005059525 


0.67032004603507 


0.67032004603564 


0.5 


0.60653065971263 


0.60653063311834 


0.60653065970940 


0.60653065971263 


0.6 


0.54881163609403 


0.54881163079527 


0.54881163609435 


0.54881163609402 


0.7 


0.49658530379141 


0.49658529860747 


0.49658530379175 


0.49658530379140 


0.8 


0.44932896411722 


0.44932896475352 


0.44932896411673 


0.44932896411722 


0.9 


0.40656965974060 


0.40656966036876 


0.40656965974011 


0.40656965974059 


1.0 


0.36787944117144 


0.36787945775656 


0.36787944117369 


0.36787944117143 



Table 2: Analytic and hybrid values of X2{t). 



t 


Analytic 


Hybrid (Af=4,M=5) 


Hybrid (Af=4,M=7) 


Hybrid (iV=4,M=9) 


0.1 


2.71451225410788 


2.71451223459708 


2.71451225410991 


2.71451225410788 


0.2 


2.45619225923395 


2.45619224196066 


2.45619225923618 


2.45619225923394 


0.3 


2.22245466204515 


2.22245467555371 


2.22245466204343 


2.22245466204515 


0.4 


2.01096013810692 


2.01096015405462 


2.01096013810538 


2.01096013810692 


0.5 


1.81959197913790 


1.81959190225067 


1.81959197912838 


1.81959197913790 


0.6 


1.64643490828208 


1.64643489597035 


1.64643490828328 


1.64643490828208 


0.7 


1.48975591137423 


1.48975590012432 


1.48975591137553 


1.48975591137423 


0.8 


1.34798689235166 


1.34798689933842 


1.34798689235054 


1.34798689235166 


0.9 


1.21970897922180 


1.21970898705545 


1.21970897922074 


1.21970897922179 


1.0 


1.10363832351433 


1.10363838023956 


1.10363832352153 


1.10363832351432 
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